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Das Sarma cond-mat/0206203| is based on a trivial mathematical mistake they have committed. 



Coulomb interaction in one-dimensional electron sys- 
tems arouses firm and enduring interest for two reasons. 
First, this is a strong interaction of clear and fundamen- 
tal nature. Second, the Coulomb interaction case is usu- 
ally difficult to treat mathematically, and therefore exact 
solutions are always important. 

We have investigated the Luttinger liquid (LL) with 
Coulomb electron-electron interaction in our paper jlj. 
The main result of our work is apparently the existence 
of a soft collective charge mode, related to dynamic elec- 
tron correlations of the 2fcp scale. Ref. Q presents also 
an analytic method to study various dynamic correlation 
functions near the threshold. 

Wang, Millis, and Das Sarma (WMS henceforth) claim 
in their Comment Q that our method, as well as ensu- 
ing results, is incorrect. The claim is based on WMS's 
abortive attempt to reproduce our calculation of the 
CDW structure factor. WMS substitute our expression 
for the structure factor into the integral equation, which 
we derived for it, and conclude that corrections to our 
result are given by diverging integrals. After that, WMS 
accuse our method to be 'inconsistent with mathematical 
analysis'. 

In this Reply, we show that WMS have committed a 
trivial mathematical mistake in their calculations. Di- 
verging integrals appear in Ref. Q just because WMS 
have incorrectly differentiated the structure factor. 

Recall that the CDW structure factor S(q, uj) contains 
the Heaviside function 9(ui — uiA as a factor, which reflects 
the existence of the threshold B We stress that 9(uj — uj q ) 
depends on uj and q, and hence must be differentiated 
when finding the derivative of S(q,uj). WMS disregard 
this fact and differentiate S(q,uj), ignoring 9- function. 
Below we demonstrate that if one acts correctly, then 
all the integrals of our paper J|] are well-defined, and 
expansions are convergent. 

In Ref. jl[] we have shown that S(q,u) satisfies the 
following integral equation: 
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Here uj q is the energy of the LL bosons, ui q — 
v F\q\/g(q), vf is the Fermi velocity, g(q) is the in- 
teraction parameter .□ For Coulomb interaction g(q) — 
/3| In l^ldl^ 1 / 2 , with d being the diameter of a quantum 
wire, (3 — [nhvp/^e 2 ] 1 ^ 2 . For simplicity, the wave num- 
ber q is measured from 2k p. 

Expand S(Q, uj — uJ q +o) in the powers of the wave num- 
ber Q, appearing in the second (frequency) argument, to 



get 
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+ ^S qq (Q,uj~uj q ) + ...}. 

Restricting expansion (J3) to the first term on the RHS, 
we have found in Ref. [flf that 
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where we denoted e = uj — u) q . The structure factor 
S(q,u>) is seen to be zero for e < 0, and to diverge as 

In our paper |l[ we have emphasized that the Coulomb 
interaction case is special in the respect that retaining 
only the first term on the RHS of Eq. (||) already gives 
the correct asymptotic behavior of S(q,ui) at e — > 0. 

On the contrary, WMS find in Ref. ^ that the second 
term on the RHS of Eq. (||) is infinite, and conclude that 
the expansion (^) does not exist at all. 

Let us prove that the second term of the expansion (||) 
is well-defined, and corrections to our result (j^) are in- 
deed small. Substitute S(q,uj) from Eq. (||) into Eq. (||), 
as WMS wish ||), but do it correctly at this time. The 
second term of the expansion becomes 
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where we denoted 6 = e — ujq. It is important that both 
e and 5 depend on q, which must be taken into account 
when calculating the derivative dan. When e < 0, 6 is 
always negative, and expression (|j) is zero due to the 
factor 0(6). If e > 0, then 6 varies from e to 0. Take 6 as 
a new integration variable. Eq. (|J) can be written as 
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As WMS note |Q| , the most important contribution to the 
integral comes from the region 6 ~ 0, where the expres- 
sion in the square brackets of Eq. (0) diverges. For this 
reason we can replace all (e — 6) in the integral with e. 
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Since the expression in the square brackets equals identi- 
cally (20)- 1 d s (exp(-4/3\lnS\ 1 / 2 )), the total formula §) 
is 
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We underline that this expression is a well-defined dis- 
tribution, rather than a classic function. Distributions 
normally arise when expanding a function with a thresh- 
old into Taylor's series. Hence it is not surprising that 
the higher order terms of the expansion diverge. Simply 
they are not 'corrections' as WMS suppose. In calculat- 
ing corrections to our result, one cannot retain only two 
terms of the expansion (g). The correct procedure |^| 
requires summing up the total series in Eq. (||) to get the 
result in terms of classic functions. We have performed 
such calculation, and have not included the result into 
our paper jlj], since it consists in a not too important 
replacement of the argument in the expression (p|) for 



S(q,iv). The argument e 
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oj — Lo q should be shifted by 
u>' q e/\ In e^/ 2 , which shift is obviously negligible as e — * 0. 

Thus we confirm our result that Eq. (^) gives the 
asymptotic behavior of S(q,ui) as e -t +0. WMS's 



claim Q that the true diverging behavior should not be 
like Eq. (||) is misleading and incorrect. 

Now consider the mistake, which is the basis of the 
WMS Comment g]. In estimating the integral (§, WMS 
first calculate the derivative of the integrand w.r.t. q. 
This way of calculation, though less economical than the 
one presented above, would nevertheless lead them to 
correct result, provided that WMS calculate the deriva- 
tive correctly. They do not take into account that the 
argument <5 = u — uj q — loq of the #(<5)-function depends 
on q, using implicitly the incorrect relation 
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which finally leads them to conclusion that the inte- 
gral (Q) diverges (see the showy, yet erroneous Eq. (14) 
of Ref. Q). In regards to WMS's appeal for the purity 
of the mathematical analysis, there is nothing left but to 
refer them to the important work of Leibniz Jd] , which ex- 
plains how the derivative of the function product should 
be found. 
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